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Loop quantum cosmology of Bianchi type II models 
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The "improved dynamics" of loop quantum cosmology is extended to include the 
Bianchi type II model. Because these space-times admit both anisotropics and non- 
zero spatial curvature, certain technical difficulties arise over and above those encoun- 
tered in the analysis of the (anisotropic but spatially flat) Bianchi type I space-times, 



and of the (spatially curved but isotropic) k=±l models. We address these and show 
that the big-bang singularity is resolved in the same precise sense as in the recent 



analysis of the Bianchi I model. Bianchi II space-times are of special interest to 
quantum cosmology because of the expected behavior of the gravitational field near 
generic space-like singularities in classical general relativity. 

PACS numbers: 98.80Qc,04.60Pp,04.60.-m 



I. INTRODUCTION 



In this paper, we will study the loop quantum cosmology (LQC) [l|, @] of the Bianchi type 
II model. These models are of special interest to the issue of singularity resolution because of 
the intuition derived from the body of results related to the Belinksii, Khalatnikov, Lifshitz 
(BKL) conjecture 0, H| on the nature of generic, spacelike singularities in general relativity 
(see, e.g., [5[]). Specifically, as the system enters the Planck regime, dynamics at any fixed 
spatial point is expected to be well described by the Bianchi I evolution. However, there are 
transitions in which the parameters characterizing the specific Bianchi I space-time change 
and the dynamics of these transitions mimics the Bianchi II time evolution. In a recent 
paper jfj], we studied the Bianchi I model in the context of LQC. In this paper we will 
extend that analysis to the Bianchi II model. We will follow the same general approach and 
use the same notation, emphasizing only those points at which the present analysis differs 
from that of [fj]. 

Bianchi I and II models are special cases of type A Bianchi models which were analyzed 
already in the early days of LQC (see in particular 0, [1]). However, as is often the case 
with pioneering early works, these papers overlooked some important conceptual and tech- 
nical issues. At the classical level, difficulties faced by the Hamiltonian (and Lagrangian) 
frameworks in non-compact, homogeneous space-times went unnoticed. In these cases, to 
avoid infinities, it is necessary to introduce an elementary cell and restrict all integrals to it 
@, 13] • The Hamiltonian frameworks in the early works did not carry out this step. Rather, 



they were constructed simply by dropping an infinite volume integral (a procedure that 
introduces subtle inconsistencies). In the quantum theory, the kinematical quantum states 
were assumed to be periodic — rather than almost-periodic — in the connection, and the 
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quantum Hamiltonian constraint was constructed using a "pre-/i " scheme. Developments 
over the intervening years have shown that these strategies have severe limitations (see, 



e.g., [Hi, ll2j, u3> ll5|)- I n this paper, they will be overcome using ideas and techniques 



that have been introduced in the isotropic and Bianchi I models in these intervening years. 
Thus, as in [6j the classical Hamiltonian framework will be based on a fiducial cell, quantum 
kinematics will be constructed using almost periodic functions of connections and quantum 
dynamics will use the "/2 scheme." Nonetheless, the space-time description of Bianchi II 
models in @, II], tailored to LQC, will provide the point of departure of our analysis. 

New elements required in this extension from the Bianchi I model can be summarized 
as follows. Recall first that the spatially homogeneous slices M in Bianchi models are 
isomorphic to 3-dimensional group manifolds. The Bianchi I group is the 3-dimensional 
group of translations. Hence the the three Killing vectors °£f on M — the left invariant vector 
fields on the group manifold — commute and coincide with the right invariant vector fields °e" 
which constitute the fiducial orthonormal triads on M. In LQC one mimics the strategy used 
in LQG and spin foams and defines the curvature operator in terms of holonomies around 
plaquettes whose edges are tangential to these vector fields. The Bianchi II group, on the 
other hand, is generated by the two translations and the rotation on a null 2-plane. Now the 
Killing vectors °£f no longer commute and neither do the fiducial triads °e\. Therefore we 
have to follow another strategy to build the elementary p la q uettes. However, this situation 
was already encountered in the k=l, isotropic models ly, |l7|. There, the desired plaquettes 



can be obtained by alternating between the integral curves of right and left invariant vector 
fields which do commute. However, in the isotropic case, the gravitational connection is 
given by A l a = c °uj l a) where °uo l a are the covectors dual to °e°: and the holonomies around 
these plaquettes turned out to be almost periodic functions of the connection component 
c [l6|, 13]. By contrast, in the Bianchi II model we have three connection components 



c 1 because of the presence of anisotropies, and, unfortunately, the holonomies around our 
plaquettes are no longer almost periodic functions of c\ (This is also the case in more 
complicated Bianchi models.) Since the standard kinematical Hilbert space of LQC consists 
of almost periodic functions of c\ these holonomy operators are not well-defined on this 
Hilbert space. Thus, the strategy |10j used so far in LQC to define the curvature operator 
is no longer viable. 

One could simply enlarge the kinematical Hilbert space to accommodate the new holon- 
omy functions of connections. But then the problem quickly becomes as complicated as full 
LQG. To solve the problem within the standard, symmetry reduced kinematical framework 
of LQC, one needs to generalize the strategy to define the curvature operator. Of course, 
the generalization must be such that, when applied to all previous models, it is compatible 
with the procedure of computing holonomies around suitable plaquettes used there. We will 
carry out this task by suitably modifying ideas that have already appeared in the literature. 
This generalization will enable one to incorporate all class A Bianchi models in the LQC 
framework. 

Once this step is taken, one can readily construct the quantum Hamiltonian constraint 
and the physical Hilbert space, following steps that were introduced in the analysis |6| of 
the Bianchi I model. However, because Bianchi II space-times have spatial curvature, the 
spin connection compatible with the orthonormal triad is now non-trivial. It leads to two 
new terms in the Hamiltonian constraint that did not appear in the Bianchi I Hamiltonian. 
We will analyze these new terms in some detail. In spite of these differences, the big bang 
singularity is resolved in the same precise sense as in the Bianchi I model [6| : If a quantum 
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state is initially supported only on classically non-singular configurations, it continues to be 
supported on non-singular configurations throughout its evolution. 

The paper is organized as follows. Section HT1 summarizes the classical Hamiltonian theory 
describing Bianchi II models. Section II III discusses the quantum theory. We first define a 
non-local connection operator A\ and use it to obtain the Hamiltonian constraint. We then 
show that the singularity is resolved and the Bianchi I quantum dynamics is recovered in the 
appropriate limit. In Section HVl we introduce effective equations for the model (with the 
same caveats as in the Bianchi I case 0]). Finally, in section V we summarize our results and 
discuss the new elements that appear in the Bianchi II model. In Appendix A we improve 
on the discussion of discrete symmetries presented in 0. The results on the Bianchi I model 
obtained in [6| carry over without any change. But the change of viewpoint is important to 
the LQC treatment of the Bianchi II model and more general situations. 



II. CLASSICAL THEORY 

This section is divided into two parts. In the first we recall the structure of Bianchi II 
space-times and in the second we summarize the phase space formulation, adapted to LQC. 



A. Diagonal Bianchi II Space-times 

Because the issue of discrete symmetries is subtle in background independent contexts, 
and because it plays a conceptually important role in the quantum theory of Bianchi II 
models, we will begin with a brief summary of how various fields are defined [l~8l . ig| . This 



stream-lined discussion brings out the assumptions which are often only implicit, making 
the discussion of discrete symmetries clearer. 

In the Hamiltonian framework underlying loop quantum gravity (LQG), one fixes an 
oriented 3-manifold M and a 3-dimensional 'internal' vector space / equipped with a positive 
definite metric qv, . The internal indices i, j, k, . . . are then freely lowered and raised by q^ 
and its inverse. A spatial triad e" is an isomorphism from / to tangent space at each point of 
M which associates a vector field v a := e*v % on M to each vector v 1 in J. 1 The dual co-triads 
are denoted by io\. Given a triad, we acquire a positive definite metric q a b : = qijUi l a ui 3 h on M. 
The metric q a b in turn singles out a 3-form e a bc on M which has positive orientation and 
satisfies e a b c edef q ad q be q c f = 3!. One can then define a 3-form on I via = e a b c e^e^el. 
Note that is automatically compatible with g^-, i.e., €ijk€i mn q d q :jm q kn = 3!. If a triad e" 
is obtained by flipping an odd number of the vectors in the triad ef, then e" and e" have 
opposite orientations and the fields they define satisfy q a b = q a b, £ a bc = tabc but = — e^. 
Had we fixed once and for all on J, then e ahc would have flipped sign under this operation 
and volume integrals on M computed with the unbarred and barred triads would have had 
opposite signs. With our conventions, these volume integrals will not change and the parity 
flips will be symmetries of the symplectic structure and the Hamiltonian constraint. 

The triad also determines an unique spin connection T l a via 

D [a ul\ = d^ + e'jkTl^ = 0. (2.1) 



1 Thus, in LQG one begins with non-degenerate triads and metrics, passes to the Hamiltonian framework 
and then, at the end, extends the framework to allow degenerate geometries. 
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The gravitational configuration variable A\ is then given by A\ = T l a + where K a ^ : = 
K l a uJbi is the extrinsic curvature of M and 7 is the Barbero-Immirzi parameter, representing 
a quantization ambiguity. (The numerical value of 7 is fixed by the black hole entropy 
calculation.) The momenta Ef carry, as usual, density weight 1 and are given by: Ef = 
v/ge". The fundamental Poisson bracket is: 

K(x), E]{ y )} = 8ttG 7 6 b J]6 3 (x,y) . (2.2) 
In Bianchi models 20, 2l|, 22|, one restricts oneself to those phase space variables admit- 



ting a 3-dimensional group of symmetries which act simply and transitively on M. Thus, 
the 3-metrics q a t, under consideration admit a 3-parameter group of isometries and M is 
diffeomorphic to a 3-dimensional Lie group G. (However, there is no canonical diffeomor- 
phism, so that there is no preferred point on M corresponding to the identity element of 
G.) To avoid a proliferation of spaces and types of indices, it is convenient to identify the 
internal space / and the Lie-algebra LG of G via a fixed isomorphism. Then, there is a 
natural isomorphism °£? between LG = I and Killing vector fields on M: for each internal 
vector v l , °^v l is a Killing field on M. For brevity we will refer to as (left invariant) 
vector fields on M. There is a canonical triad °e" — the right invariant vector fields — which 
is Lie dragged by the This triad and the dual co-triad °uf a satisfy: 

[%,%} = 0, [°e l ,%} = -°C* e k , 
L„ Ci (V) =0, d °u k = \ °CJV A V, (2.3) 



where °C% denotes the structure constants of LG. It is convenient to use the fixed fields °ef 
and °uj l a as fiducial triads and co-triads. 

In the case when G is the Bianchi II group, we have °C\ k = as in all class A Bianchi 
models and, furthermore, the symmetric tensor k kl := e^ 1 has signature +,0,0. There- 
fore, we can fix, once and for all an orthonormal basis °b\ 1 b\ 1 °b\ in I such that the only 
non-zero components of °G\, are 

"Cas = -CI2 = « , (2-4) 
where a is a non-zero real number. 2 We will assume that this basis is so oriented that 

em := £ijk°b\%% = e (2.5) 

where e = ±1 depending on whether the frame e a { (which determines the sign of e^) is right 
or left handed. Throughout this paper we will set = °^1°b\, °e\ = °e1°b\, %j\ = °uj % a °b\, 
etc. 

The form of the components of in this basis implies that M admits global coordinates 
x, y, z such that the Bianchi II Killing vectors have the fixed form 



2 Without loss of generality a can be chosen to be 1. We keep it general because we will rescale it later 
(see Eq. (|2.17[) ) and because we want to be able to pass to the Bianchi I case by taking the limit a — >• 0. 
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These expressions bring out the fact that, if we were to attempt to compactify the spatial 
slices to pass to a T 3 topology — as one can in the Bianchi I model — we will no longer have 
globally well-defined Killing fields. Thus, in the Bianchi II model, we are forced to deal with 
the subtleties associated with non-compactness of the spatially homogeneous slices. 
In the x, y, z chart, the right invariant triad is given by 



*- fc (sH£)'- *-(!)'■ (2 - 7 » 

and the dual co-triad by 

%)\ = {dx) a - &z(dy) a , °u 2 a = (dy) a , lu 3 a = (dz) a . (2.8) 

They determine a fiducial 3- metric °q a b '■— qifulfal with Bianchi II symmetries: 

°q ab dx a dx b = (dx - az dy) 2 + dy 2 + dz 2 . (2.9) 

In the diagonal models, the physical triads e" are related to the fiducial ones by 3 

c^ = a*(T)K, and o^e? = X (2-10) 

where the a, are the three directional scale factors. Since the physical spatial metric is given 
by q a b = u l a oj bi) the space-time metric can be expressed as 

ds 2 = -Ndr 2 + ai(r) 2 (dx - az dy) 2 + a 2 (r) 2 dy 2 + a 3 (r) 2 d^ 2 (2.11) 

where N is the lapse function adapted to the time coordinate r. 

For later use, let us calculate the spin connection (12. ip determined by triads ef. From 
the definition of it follows that 

r*« = ej U a uj b]k + l -etu l a d [c u h] ^ . (2.12) 

Using (I2.5p . the components of Y l a in the internal basis °b\, °b l 2) °b l 3 can be expressed as 

ri = -^Uj ; r 2 = -^i^ 2 ; r 3 = -^i^ 3 . (2.13) 

a 2a 2 a 3 a a 2a 3 2a 2 a V ' 

Before studying the dynamics of the model, let us examine the action of internal parity 
transformation 11^ which flips the kth triad vector and leaves the orthogonal vectors alone. 
(For details see Appendix and [19]). Under the parity transformation n 1? for example, we 
have: e\ — > — e", e 2 — > e 2 , e\ — > e 3 and a% — > —a\, a 2 — > a 2 ,a 3 — ► a 3 whence — > 
— T", T 2 — ► T 2 , T 3 — > r|. Thus, both ef and are proper internal vectors, e on the other 
hand is a pseudo internal scalar, e — > — e under every 11^. Note that the fiducial quantities 
carrying a label o do not change under this transformation; it affects only the physical 
quantities. 



3 There is no sum if repeated indices are both covariant or contravariant. As usual, the Einstein summation 
convention holds if a covariant index is contracted with a contravariant index. 
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B. The Bianchi II Phase space 

As is usual in LQC, we will now use the fiducial triads and co-triads to introduce a 
convenient parametrization of the phase space variables, Ef, A\. Because we have restricted 
ourselves to the diagonal model and these fields are symmetric under the Bianchi II group, 
from each equivalence class of gauge related phase space variables we can choose a pair of 
the form 

E? = PiVWVi and K = (2-14) 

where, as spelled out in footnote 3, there is no sum over i. Thus, a point in the phase space 
is now coordinatized by six real numbers pi,c l . One would now like to use the symplectic 
structure in full general relativity to induce a symplectic structure on our six- dimensional 
phase space. However, because of spatial homogeneity and the M 3 spatial topology, the 
integrals defining the symplectic structure, the Hamiltonian (and the action) all diverge. 
Therefore we have to introduce a fiducial cell V and restrict integrals to it 

BOH. We will 

take the fiducial cell to be rectangular with edges along the coordinate axes and lengths of 
Li, L 2 and L 3 with respect to the fiducial metric °q a b. It then follows that the volume of the 
fiducial cell with respect to °q a b is V a = L\L 2 L 3 . Then the non-zero Poisson brackets are 
given by: 

where 7 is the Barbero-Immirzi parameter. As in the Bianchi I case, we have a 1-parameter 
ambiguity in the symplectic structure because of the explicit dependence on V Q and we 
have to make sure that the final physical results are either independent of V Q or remain 
well-defined as we remove the 'regulator' and take the limit V — > 00. 

It is convenient to rescale variables to absorb this dependence in the phase space coordi- 
nates fas was done in the treatment of Bianchi I model in 61). Let us set 



{c\P 3 } = (2.15) 



Pi = L 2 L 3 pi, p 2 = L 3 L ± p 2 , p 3 = L ± L 2 p 3 , (2.16) 

ci = L x ci, c 2 = L 2 c 2 , c 3 = L3C3 and a = ^ 2 ^ z a , (2.17) 

where the last rescaling has been introduced to absorb factors of L, which would otherwise 
unnecessarily obscure the expression of the Hamiltonian constraint. The Poisson brackets 
between these new phase space coordinates is given by 

{c\ Pj } = ZnG^d). (2.18) 

These variables have direct physical interpretation. For example, pi is the (oriented) area of 
the 2-3 face of the elementary cell with respect to the physical metric q a b and = exp c\T\ 
is the holonomy of the physical connection A l a along the first edge of the elementary cell. 

Our choice (12.141) of physical triads and connections has fixed the internal gauge as well 
as the diffeomorphism freedom. Furthermore, it is easy to explicitly verify that, thanks to 
f)2.14p . the Gauss and the diffeomorphism constraints are automatically satisfied. Thus, as 
in[6|, we are left just with the Hamiltonian constraint 

C H = f \ ^ E ' E U e^ k F ab k - 2(1 + l 2 )K\ a K^ + A^ matt ] d 3 x , (2.19) 
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F ab k = 2d [a A k b] + e io k A\A{ (2.20) 



where 



is the curvature of A l a and 7i ma tt is the matter Hamiltonian density. As in [6j, our matter 
field will consist only of a massless scalar field T which will later serve as a relational time 
variable a la Liebniz. (Additional matter fields can be incorporated in a straightforward 
manner, modulo possible intricacies of essential self-adjointness.) Thus, 

T^matt = 7, ^J-t - (2-21) 



Since we want to use the massless scalar field as relational time, it is convenient to use a 
harmonic-time gauge, i.e., assume that the time coordinate r in (12.111) satisfies Or = 0. The 
corresponding lapse function is iV = \Z\pijhP3~\- With this choice, the Hamiltonian constraint 
simplifies considerably. Note first that the basic canonical variables can be expanded as 

E« = ^U^\X and K = jjy» (2.22) 
and the extrinsic curvature is given by 

K = r\A\-r a ). 

Next, using p\ = (sgnai) |a2«3 1 L2L3 etc, the components of the spin connection become: 

2pj Li 2pi L 2 2pi L 3 

Collecting terms, the Hamiltonian constraint (12.191) becomes 

1 



C 



H 



8nGr 



P1P2C1C2 + P2P3C 2 c 3 + P3P1C3C1 + aep 2 p3Ci 



-(l + 7 2 )(^?) 2 +^ (2.24) 



2pi 



1 



2 



r m 1 



' H 87rG7 2 



aep 2 p 3 ci -I + 7 (— J , 2.25 

2pi -I 



where Cjj is the Hamiltonian constraint (including the matter term) for Bianchi I space- 
times which has already been studied in [6f. Note that this constraint is recovered in the 
limit a — > 0, as it must be. 

Knowing the form of the Hamiltonian constraint, it is now possible to derive the time 
evolution of any classical observable O by taking its Poisson bracket with Ch'- 

6 = {<D,Ch}, (2.26) 

where the 'dot' stands for derivative with respect to harmonic time r. This gives 

Pi = l~ l {p\Pic<i + P1P3C3 + aep 2 p 3 ), (2.27) 

P2 = 1~ 1 {P2P\C\ +P2P3C3), (2.28) 
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Pa = 1 1 {P3PiCi +P3P2C2), (2.29) 
ci = -i Uc lC2 + P3Cl c 3 + -L(l + 7 2 ) (^) 2 ) , (2.30) 

C2 = Pic 2 ci +P3C2C3 + «£p 3 ci - — (1 + 7 )( ) ), (2-31) 

7 V 2p 2 Pi ' 

c 3 = — fpit&ci +P2C3C2 + aep 2 ci - -^-(1 + 7 2 ) ( ap2P3 ) 2N ) . (2.32) 
7 V 2p 3 p x ' J 

Any initial data satisfying the Hamiltonian constraint can be evolved by using the six equa- 
tions above. It is straightforward to extend these results if there are additional matter 
fields. 

Finally, let us consider the parity transformation 11^ which flips the /cth physical triad 
vector e a k . (As noted before, this transformation does not act on any of the fiducial quantities 
which carry a label o.) Under this map, we have: q a b — > q a b, tabc ^abc but — > e — ► 

— e. The canonical variables c l ,pi transform as proper internal vectors and co-vectors: For 
example 

n 1 (ci,c 2 ,c 3 ) -> (-ci,c 2 ,c 3 ) and IIi(pi,p2,p 3 ) ->• (-Pi,P2,p 3 ) ■ (2.33) 

Consequently, both the symplectic structure and the Hamiltonian constraint are left invari- 
ant under any of the parity maps 11^. 

This Hamiltonian description will serve as the point of departure for loop quantization 
in the next section. 

III. QUANTUM THEORY 

This section is divided into three parts. In the first, we discuss the kinematics of the 
model, in the second we define an operator corresponding to the connection A\ using 
holonomies and in the third we introduce the Hamiltonian constraint operator and describe 
its action on states. 

A. LQC Kinematics 

The kinematics for the LQC of Bianchi II models is almost identical to that for Bianchi 
I models. Therefore, in the sub-section we closely follow [6]. 

Let us begin by specifying the elementary functions on the classical phase space which 
will have unambiguous analogs in the quantum theory. As in the Bianchi I model, the 
elementary variables are the momenta pi and holonomies of the gravitational connection A\ 
along the integral curves of the right invariant vector fields °e°\ Let Tj be a basis of the Lie 
algebra of SU(2), satisfying TjTj = \tij k Tk — t<%I where I is the unit 2x2 matrix. Consider 
an edge of length £L k with respect to the fiducial metric °q a b, parallel to °e^. The holonomy 
hjjjr along it is given by 

h ( j:\ci, c 2 , c 3 ) = exp (£c k r k ) = cos ^1 + 2 sin ^r k . (3.1) 
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(Note that £ depends of the fiducial cell but not on the fiducial metric.) This family of 
holonomies is completely determined by the almost periodic functions exp(i£ck) of the con- 
nection. These almost periodic functions will be our elementary configuration variables 
which will be promoted unambiguously to operators in the quantum theory. 

It is simplest to use the p-representation to specify the gravitational sector 7if™ V °f the 
kinematic Hilbert space. The orthonormal basis states \pi1P21P3) are eigenstates of quantum 
geometry. For example, in the state \pi,P2,P3) the face S23 of the fiducial cell V (given by 
x =const) has area \p\\. The basis is orthonormal in the sense 

(Pl,P2,P3\p'l,P 2 ,P3) = ^ PlP '^p 2P '^p 3 p' 3 , (3-2) 

where the right side features Kronecker symbols rather than the Dirac delta distributions. 
Hence kinematical states can consist only of countable linear combinations 

I*) = *(fi>P2,P3)bi,P2,P3) (3.3) 

P1,P2,P3 

of these basis states for which the norm 

\\m 2 = \^(puP2,P3)\ 2 (3.4) 

P1,P2,P3 

is finite. Because the right side features a sum over a countable number of points on 
IR 3 , rather than a 3-dimensional integral, LQC kinematics are inequivalent to those of the 
Schrodinger approach used in Wheeler-DeWitt quantum cosmology. 

Next, recall that on the classical phase space the three reflections IT : e" — > —el are 
large gauge transformations under which physics does not change (since both the metric and 
the extrinsic curvature are left invariant). These large gauge transformations have a natural 
induced action, denoted by IT, on the space of wave functions \E f (p 1 ,p 2 ,P3)- For example, 

ni*(pi,p2,p 3 ) = ^(-p 1 ,p2,p 3 )- (3.5) 

Since TI? is the identity, for each i, the group of these large gauge transformations is simply 
Z2. As in Yang- Mills theory, physical states belong to its irreducible representation. For 
definiteness, as in the isotropic and Bianchi I models, we will work with the symmetric 
representation. It then follows that is spanned by wave functions ^(pi,P2,P3) which 

satisfy 

*(Pl,P2,P3) =*(blUP2|,bs|) (3-6) 

and have a finite norm (13.41) . 

The action of the elementary operators on Ti.f^ is as follows: the momenta act by mul- 
tiplication whereas the almost periodic functions in shift the ith argument. For example, 

[Pi^]{pi,P2,P3) = Pi *(pi,P2,Ps) and exp(i£ci)^ {pi,p 2 ,P3) = *(pi - 87i-fGH£,p 2 ,P3) ■ 

_ (3.7) 

The expressions for p 2 , exp(i£c 2 ),P3 and exp(i£c 3 ) are analogous. Finally, we need to define 
the operator e since e features in the expression of the Hamiltonian constraint. In the 
classical theory, e is unambiguously defined on non-degenerate triads, i.e., when P1P2P3 7^ 0. 
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In quantum theory, wave functions can have support also on degenerate configurations. We 
will extend the definition to degenerate triads using the basis |pi,P2,P3) : 

\pi,P2,Pa) i f P1P2P3 > 0, 
£\Pi,P2,Pa) ■= i . . n (3.8) 

" |Pl,P2,P3) II PlP2P3 < 0. 



grav 
kin 



Finally, the full kinematical Hilbert space H^m will be the tensor product H^ n = H 

where Wjgf* = L 2 (1R, dT) is the matter kinematical Hilbert space for the homogeneous 
scalar field. On 7^™„ tt , T will act by multiplication and pr '■= —ihdr will act by differentia- 
tion. As in isotropic and Bianchi I models, our final results would remain unaffected if we 
use a "polymer representation" also for the scalar field. 



B. The connection operator A l a 

To define the quantum Hamiltonian constraint, we cannot directly use the symmetry 
reduced classical constraint (12.241) because it contains connection components themselves 
and in LQC only almost periodic functions of Ch have unambiguous operator analogs. Indeed, 
in all LQC models considered so far 0, [Io|, 0, H, E3, 23, 0, 25[, we were led to return to 



the expression f 1 2 . 1 9 j) in the full theory and mimic the procedure used in LQG [26|]. More 
precisely, the key strategy was to follow full LQG (and spin foams) and define a "field 
strength operator" using holonomies around suitable closed loops. In the Bianchi I model, 
these closed loops were formed by following integral curves of right invariant vector fields 
(which are also left invariant). As mentioned in section [Til i n the Bianchi II model the 
right invariant vector fields define the fiducial triads °e", the left invariant vector fields, the 
Killing fields °£\ Neither constitutes a commuting set, whence their integral curves cannot 
be used to form closed loops. However, as in the k=l case Il7l |. one can hope to exploit 



the fact that the right invariant vector fields do commute with the left invariant ones and 
construct the closed loops by alternately following right and left invariant vector fields. But, 
as mentioned in section [Tj a new problem arises: unlike in the k=l (or Bianchi I) model the 
resulting holonomies are no longer almost periodic functions of Ck, whence the Hilbert space 
7ikin V does not support these holonomy operators. For completeness we will first show this 
fact explicitly and then introduce a new avenue to bypass this difficulty. 

The problematic curvature component turns out to be F yz l . To^ construct the corre- 



sponding operator, following the strategy used in the k=l case [ly, |17||, we will construct a 
closed loop D yz as follows. In the coordinates (x,y,z), i) Move from (0,0,0) to (0,p, 2 L 2 ,0) 
following £2; ii) then move from (0,fl 2 L 2 ,0) to (0,p, 2 L 2 ,p, 3 L3) following °e\\ iii) then move 
from (0, fi 2 L 2 , /X3L3) to (0, 0, //3L3) following — £f; and, finally, iv) move from (0,0, //3L3) 
to (0, 0, 0) following —°e%. The parameters /2, which determine the 'lengths' of these edges 
can be fixed by the semi-heuristic correspondence between LQC and LQG exactly as in the 
Bianchi I model 0] because the geometric considerations used in that analysis continue to 
hold without any modification in all Bianchi models with M 3 spatial topology: 



(1| NH, / jgsjAgi 

Vl IP2P3 1 ' 2 V \pips\ ' 3 V \pm\ 

where A£p { = 4-^/3^7 £p { is the 'area gap'. The holonomy around this closed loop D yz is 
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given by 

h D = _ 2 cos f sin (^) (c 2 sin(// 3 c 3 ) + a/2 3 ci cos(// 3 c 3 )) (3.10) 
c uoUzLoLz V 2 / V 2 / V J 



c ii 2 ^L 2 L 

where 



c= Ja^lcl + cl (3.11) 



If we were to shrink the loop so that the area it encloses goes to zero, we do indeed recover 
the classical expression of F yz l . However, because of presence of the term c, if a ^ the 
right side fails to be almost periodic in c\ and c 2 - Hence this holonomy operator fails to 
exist on 7ikin- It is clear from the expression (13.111) of c that the problem is independent of 
the specific way /2« are fixed. 

We will bypass this difficulty by mimicking another strategy used in full LQG [26]: We 
will use holonomies along segments parallel to °e a i to define an operator corresponding to 
the connection itself. This is a natural strategy because holonomies along these segments 
suffice to separate the Bianchi II connections (I2.14p . Let us set A a := A k a T k . Then we have 
the identity: 

^&£^K >, -(^>)-') (3.12) 

k 

where h£ is given by (13. ip . In the expressions of physically interesting operators such 
as the Hamiltonian constraint of full LQG, one often replaces A a with the (analog of the) 
right side of (I3.12j) . But because of the specific forms of these operators, the limit trivializes 
on diffeomorphism invariant states of LQG. In LQC, we have gauge fixed the system and 
therefore cannot appeal to diffeomorphism invariance. Indeed, while the holonomies are 
well-defined for each non-zero £ k , the limit fails to exist on Wjj^. A natural strategy is to 
shrink l k to a judiciously chosen non-zero value. But what would this value be? In the 
case of plaquettes, we could use the interplay between LQG and LQC directly because the 
argument Pi of LQC quantum states refers to quantum areas of faces of the elementary cell V 
0]. For edges we do not have such direct guidance. There is, nonetheless a natural principle 
one can adopt: Normalize l k such that the numerical coefficient in front of the curvature 
operator constructed from the resulting connection agrees with that in the expression of the 
curvature operator constructed from holonomies around closed loops, in all cases where the 
second construction is available. We will use this strategy. Let us apply it to the Bianchi I 
model where F a b h = €ij k A % a A\. Using holonomies around closed loops one obtains the field 
strength operator 

Fj = e l3 k y a y UY (3.13) 

where 

, sin fjc Q y xsin/^q ^ , 
{-—— LU a ) = {— - — u a ) (no sum over l) 
\lL Hi-Li 



see Eqs (3.12) and (3.13) in [6J). Therefore, our strategy yields £k = 2/z fc , that is, 



% = ^j^lt, (3-14) 



sin(/2 fc c fc ^ 

where there is no sum over k. Note that the principle stated above leads us unambiguously 
to the factor 2 in £ k = 2fl k ; without recourse to a systematic strategy, one may have naively 
set 4 = p k . 
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If we compare the expression (13.141) of the connection operator with the expression (I2.14p 
of the classical connection, we have effectively defined an operator c via 



Cfe 



sm([i k c k ) 



(3.15) 



where there is again no sum over k. In the literature such a quantization of c is often called 
"polymerization." Our approach is an improvement over such strategies in two respects. 
First, we did not just make the substitution c — > sin£c/£ by hand; a priori one could have 
used another substitution such as c — > tan£c/£. Rather, as in full LQG, we used the strategy 
of expressing the connection in term of holonomies, 'the elementary variables'. But this still 
leaves open the question of what £ one should use. We determined this by requiring that 
the overall normalization of F a b k constructed from A l a = d {L' l )~ 1 °uj l a should agree with 
that of F a b k constructed from holonomies around appropriate closed loops, when the second 
construction is possible. Therefore, our construction is a bona-fide generalization of the 
previous constructions used successfully in LQC. 

This strategy has some applications beyond the Bianchi II model studied in this paper. 



24j . The analysis uses the [i 



First, the k=— 1 isotropic case has been studied in detail in [2c 
scheme, carries out a numerical simulation usin g ex act LQC equations and shows that the 
effective equations of the "embedding approach" [27], [28| (discussed in section HV|) provide an 
excellent approximation to the quantum evolution. While this is an essentially exhaustive 
treatment, as 23j, |24j itself points out, the treatment has a conceptual limitation in that it 
builds holonomies around the closed loops using the extrinsic curvature K l a — rather than 
A\ — as a "connection" . This limitation can be overcome in a straightforward fashion using 
our current strategy. More importantly, this strategy is applicable to all class A Bianchi 
models, including type IX. Thus, it opens the door to the LQC treatment of all these models 
in one go. 



C. The quantum Hamiltonian constraint 

With the connection operator at hand, one can construct the Hamiltonian constraint 
operator starting either from the general LQG expression (12.191) or the symmetry reduced 
expression (12.241) . We will begin by a small change in the representation of kinematical 
states which will facilitate this task. 



by 



1. A more convenient representation 
Ignoring factor ordering ambiguities for the moment, the constraint operator Ch is given 



C 



1 



H 



8%G-? 2 A£ 2 



pi 



P1P2IP3I sin/iici sin/i 2 C2 + \p1\p2P3 sin/i 2 c 2 sin /i 3 c 3 

1 



+ Pi\P2\P3 sin /i 3 c 3 sin [L X c x 
2 



87rG 7 2 



+ - 2 Pt 



(3.16) 
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where for simplicity of notation here and in what follows we have dropped the hats on the 
Pi and sin /ijCj operators. Recall that, classically, the Bianchi II symmetry group reduces to 
the Bianchi I symmetry group if we set a = 0. If one sets a = in (13 . 16j) . the last two terms 
disappear and the operator C H reduces to that of the Bianchi I model 0] showing explicitly 
that our construction is a natural generalization of the strategy used there. 

To obtain the action of operators corresponding to terms of the form sin /ZjCj we use the 
same strategy as in [6]. As shown there, it is simplest to introduce dimensionless variables 

= sgn^y^ . u) 

1 (47r 7 v / A4 1 ) 1 / 3 ' 

Then the kets I2, h) constitute an orthonormal basis in which the operators pk are diag- 
onal 

Pk\h,h,h) = [sgn(4)(47r 7 v / A4i) 2/3 ^ \h,h,h). (3.18) 

Quantum states will now be represented by functions ^(li,h,h)- The operator e^ lCl acts 
on them as follows 

[e 1 ^ *](h,l2,h) = *(h- Jn,h,h) 

1*2*3 1 
V 

where we have introduced the variable v = 2/1/2/3 which is proportional to the volume of 
the fiducial cell: 

V*(h,l 2 ,k) = [2Tr-yVE\v\4A*(h,l2,h). (3.20) 

(The e*^ ici operator is well-defined in spite of the appearance of \fah\ in the denominator; 
see Q.) The operators e^ 2C2 and e lp,3C3 have analogous action. 

We are now ready to write the Hamiltonian constraint explicitly in the /^-representation. 
As noted above, the three terms in the first square bracket on the right hand side of Eq. 
(I3.16P constitute the gravitational part of Ch for the LQC of Bianchi I model 4 and have 
been discussed in j6[. In the next two sub-sections we will now discuss the last two terms, 
which are specific to the Bianchi II model. 



2. The Fourth term in Ch 
Using a symmetric factor ordering, the fourth term becomes 

dg) = - Q f^V» \^ { e^ lCl + ^ lCl e)\^. (3.21) 



167rG7 2 \/A4>i 

(Note that P2 and ^3 commute with the other terms in ). The operator p\ is self-adjoint 
on whence any measurable function of pi is also a well-defined self-adjoint operator. 

However, since kets |/i = OJ2J3) are normalizable in H^ v , the naive inverse powers of p\ 



4 There are some minor changes in the action of these three terms since 7 is no longer treated as a pseu- 
doscalar (see Appendix EJ , but these do not affect the discussion. 
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fail to be densely defined and cannot be self-adjoint. To define inverse powers, as is usual in 
LQG, we will use a variation on the Thiemann inverse triad identities (26[. Classically, we 
have the identity 



Pi 



-1/4 



isgnbO /bgpal e -^c 1{ ^ lcsbi|1 /4 }- (322) 



27rG 7 V A£ 2 



pi 



which holds for any choice of p,\. Since it is most natural to use the same jli that featured 
in the definition of the connection operator, we will make this choice. Eq f 1 3 . 2 2 j) suggests a 
natural quantization strategy for |pi|~ 1//4 . Using it and the parity considerations, we are led 
to the following factor ordering 



5 



^ = -^f^v^ e ~ i " ici/2 k [emici > bi|1/4] e_i " ici/2 ' (3 - 23) 

where, as is common in LQC, sgn(pi) is defined as 

( +1 if pi >0, 

sgn( Pl )=J Oifp 1 = 0, (3.24) 
{ -1 if p! < 0. 

At first it may seem surprising that the expression of |pi|~ 1//4 involves operators other 
than pi. It is therefore important to verify that it has the standard desirable properties. 
First, as one would hope, it is indeed diagonal in the eigenbasis of the operators pk~- 



(« 3 )l \h,h,h)- 



(47r7v / A£| 1 ) 1 /6 V 

(3.25) 

Second, on eigenkets with large volume, the eigenvalue is indeed well-approximated by p l 1 ^ 4 , 
whence on semi-classical states it behaves as the inverse of p 1 / 4 , just as one would hope. 

Thus, (13.251) is a viable candidate for |pi|~ 1//4 . But there are interesting non-trivialities in 
the Planck regime. In particular, although counter-intuitive, as is common in LQC the 
operator annihilates states \h,h, h) with v = 2lil 2 l 3 = 

Finally, note that the operator e appearing in the expression (13.211) of 6$ either operates 

immediately before or after |pi| -1 / 4 . Since |pi| -1 / 4 annihilates all zero volume states and 
e acts on such states as the identity operator, we only need to consider the action of e on 
states with nonzero volume. In this as sgn(f ). Therefore the action of can 

be written as: 

>(4) , T /| /, , t \ iany/Aht^! , w , 4 
C * *J k ' h) = - ( 47r7 VA)Va s ^ ^ 

\v + sgn(l 2 l 3 )\ - y/\v-sgn(l 2 l 3 )\) \$ + {h, k,h) - $"(Zi, h, h)} (3.26) 



5 In the classical theory, (L2L3) 1 / 4 |pi| -1 ^ 4 is independent of the choice of the elementary cell. As pointed 
out in 23] the inverse triad operators, by contrast, depend on the choice of the cell. However, one can 
verify that as we remove the regulator, i.e., take the limit V — > K 3 , as in the classical theory, the limiting 
(L 2 i3) 1/4 IpiP^ 4 has a well defined limit. 
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where 

$Hh,h,k) =(y/\v ± 2sgn(Z 2 / 3 ) + sgn(Z 2 Z 3 ))| - y/\v ± 2sgn(/ 2 Z 3 ) - sgn(Z 2 Z 3 ; 

x (sgn(t;) + sgn(t; ± 2sgn(/ 2 / 3 ))) ; / 3 ). (3.27) 

Recall that in the classical theory the singularity corresponds precisely to the phase 
space points at which the volume vanishes. Therefore, as in the Bianchi I model, states with 
support only on points with v — will be called 'singular' and those which vanish at points 
with v — will be called regular. The total Hilbert space Ti^ r is naturally decomposed 
as a direct sum = Ti^g ® ^feg V °f singular and regular sub-spaces. We will conclude 

this discussion by examining the action of C# on these sub-spaces. Note first that in the 

action (I3.2ip of C^, the state is first acted upon by the operator |pi|" 1//4 . Since this operator 
annihilates states |/i/ 2 ,Z 3 ) with v = 2ZiZ 2 Z 3 = 0, singular states are simply annihilated by 
Cjp. In particular this implies that the singular sub-space is mapped to itself under this 
action. It is clear from (13.271) that if \1/ is regular, i.e. vanishes on all points with v = 0, 
$ ± also vanish at these points. Thus the regular sub-space is also preserved by this action. 
This fact will be used in the discussion of singularity resolution in section IIII C4[ 

Remark: Our definition of the operator |p| -1 / 4 is not unique; as is common with non- 
trivial functions of elementary variables, there are factor ordering ambiguities. For example, 
for < n < 1/2, we have the classical identity 



I r i/ 2 = -^gn(pQ^3l e ^ lcl { * } 



Hence, it is possible to instead define p l 1 ' i as 

-l/(4n-2) 



pi 1/4 = (bTp/ 2 ) 




(47r 7v / A4 1 ) (2+2n)/3 



4 n (8ir<yGVAe Pl ) 3 n 



sgn(/ 1 )|/ 2 / 3 | 1 -^ \v + sgn(/ 2 Z 3 )| 2 " - \v - sgn(Z 2 / 3 )| 2 " 



For n 7^ 1/4, this choice for the operator p l ^ 4 is not equivalent to the one we chose. 
These two choices are both well-defined and admit the same classical limit but they differ in 

the Planck regime. It is also possible to construct other such inequivalent p[ X ^ candidate 
operators. For definiteness we have made the 'simplest' choice. 

3. The fifth term in Ch 

Let us now consider the last term in the expression of the gravitational part of the 
Hamiltonian constraint 

^ ) = 32^ (1 + 72)( ^ 3)2 ^ 2 - (3 ' 28) 
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This term is simpler since it only involves powers of pk and we are working in a representation 
where pk are diagonal. From our discussion of the last section, it is natural to set 

^2 ._ ( "~-L/4 



pr 2 :=bi ) , (3-29) 

then we have 



i ,-„ iTi/; ? ? _8^A(l + 7 ^ 

(4tt 7 \/A) 2 / 3 V ; 3 



x 



[yj\v + sgn(Z 2 Z 3 )| - yj\v - sgn(Z 2 Z 3 )|) *{h,h,h). (3.30) 



Again, it is clear that if v — 0, the wave function is annihilated by this part of the constraint. 
Also, it follows by inspection that the singular and regular subspaces are both mapped to 
themselves by the action of Cjp. 



4- Singularity resolution 

We can now determine the gravitational part C grav of the Hamiltonian constraint by 
combining the results of [g[ and Eqs. (I3.26p and f)3.30p . We have: 

Cgrav = Cgfav + + (3.31) 

where Qrav is the gravitational part of the Hamiltonian constraint in the Bianchi I model 
0]. There is however, a conceptual subtlety. In the classical theory the Hamiltonian density 
Cgra.v/(LiL 2 L 3 ) 2 is independent of the choice of the elementary cell (where we have to divide 
by (L1L2Z/3) 2 because the lapse corresponding to harmonic time scales as [L\L 2 L 3 ) and the 
Hamiltonian constraint is obtained by integration over the elementary cell V). As shown in 

the section V of [6], C g fav / '(LiL 2 L 3 ) 2 is again independent of the choice of the elementary 
cell V. However, the two additional terms that are special to the Bianchi II model are not 



independent of V because they involve the inverse-triad operators 23|]. Nonetheless, in the 
limit as we take the regulator away, i.e., V — > M 3 , the operator C grav / '(LiL 2 L 3 ) 2 has a well- 
defined limit (see footnote 5). Strictly speaking, in the discussion of Bianchi II quantum 
dynamics, we have to work with this limit, rather than with operators defined using a fixed 
cell. 

As in the Bianchi I model, the action simplifies if we replace one of the k by v. In the 
Bianchi I model, it does not matter which of the U is replaced because of the additional 
symmetry of that model. In the Bianchi II case, while it remains possible to replace any of 
the li, it is simplest to replace l\ by v and represent quantum states as $ = *$?(l 2 , l 3 , v ; T). 
This change of variables would be nontrivial if, as in the Wheeler-DeWitt theory, we had used 
the Lesbegue measure in the gravitational sector. However, it is quite tame here because 
the norms are defined using a discrete measure. The inner product on T^irT ^ s now gi ven by 

(*i|^ 2 )kin=^ *i{h,h,v)Mh,h,v) (3.32) 

h,l3,V 

and states are symmetric under the action of Ilk- In Appendix |Aj we show that, under 
the action of reflections n^, the operators sin/IjCi have the same transformation properties 
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that q have under reflections Ilj in the classical theory. As a consequence, C grav is also 
reflection symmetric. Therefore, its action is well defined on 7if™ v : C grav is a densely defined, 
symmetric operator on this Hilbert space. In the isotropic case, its analog has been shown 
to be essentially self-adjoint j30|. In what follows we will assume that (13.311) is essentially 
self-adjoint on TC^'a and work with its self-adjoint extension. 

It is now straightforward to write down the full Hamiltonian constraint on Hf^: 

-h 2 d*ty(l 2 J 3 ,v;T) =QV(l 2 ,l 3 ,v;T) where = -2C grav (3.33) 



As in the isotropic case [29(, one can obtain the physical Hilbert space 7i p h y by a group 
averaging procedure and the final result is completely analogous. Elements of 7i p h y consist 
of 'positive frequency' solutions to (13.331) . i.e., solutions to 

-iHd T V(l2,h,v;T) = ^/\e\*(h,h,v;T), (3.34) 

which are symmetric under the three reflection maps Ilj, i.e. satisfy 

*(! 2 ,l 3 ,v-,T) = *(\l 2 \,\h\,H,T). (3.35) 

The scalar product is given simply by: 

(^i|^ 2 ) P hys = (^i(k,h,v;T )\^ 2 (l 2 ,l 3 ,v;T )) UD 

= MWM^To) (3.36) 

h,l2,h 

where T is any "instant" of internal time T. 

We can now address the issue of singularity resolution using general properties of various 
operators. Recall that the gravitational part of the Hamiltonian constraint operator in the 
Bianchi I model shares two properties with the fourth and the fifth terms studied above 
which are specific to the Bianchi II model. First, it annihilates singular states and, second, 
singular states decouple from the regular states under its action. Therefore the full Bianchi 
II Hamiltonian constraint also has these two properties. Since the singular states decouple 
from regular states 6 , an initial state in the regular sub-space cannot become singular during 
evolution. It is in this precise sense that the classical singularity is resolved. Sometimes 
one considers weaker forms of singularity resolution. For example, it could happen that 
the evolution of the wave function is always well defined but a regular state can evolve to 
the singular sub-space. For the Bianchi I and II models, the singularity is resolved in a 
stronger sense: Not only is the evolution well defined at all times, but the singular states 
(are stationary and) decouple entirely from the regular ones. 



5. The explicit form of the Hamiltonian constraint 

We will conclude by providing an explicit form of the full quantum constraint equation 
that will be needed in numerical simulations. 



6 Singular states are in the kernel of 8 and regular states are orthogonal to the singular ones. From spectral 
decomposition one expects \f@ to have the same property. However, to complete this argument, one would 
have to establish that C gYav is essentially self-adjoint and its self adjoint extension also shares this property. 
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Recall that in the Bianchi I model [6j symmetries enabled us to restrict our attention 
to the positive octant of the 3-dimensional space spanned by {lx,l 2 ,l 3 )- This is again the 
case for the Bianchi II model. More precisely, elements of 7^|[^ v are invariant under the 
three parity maps Ilk and, as shown in the Appendix |Aj the Hamiltonian constraint satisfies 
fifcCgravIIfc = C grav . Therefore, knowledge of the restriction of the image Cg rav ^ of $ to the 
positive octant suffices to determine C gra v^ completely. In the positive octant, sgn(/ fc ) can 
only be or 1 which simplifies the action of operators. Therefore, in the remainder of this 
section we will restrict the argument of Ch^ to the positive octant. The full action is given 
simply by 

(4rav*)(*2,M) = (4rav*)(N,M>|)- (3-37) 



Since the singular states are annihilated by C gTav , their evolution is trivial: 

d*q(l 2 ,l 3 ,v = 0;T) = Q. 



(3.38) 



Non-singular states are physically more relevant. On them, the explicit form of the full 
constraint is given by: 



tvG 



d*y{l 2 ,h,v;T)=— 



^{{v + 2)v^+4 +(Z 2 , l 3 , v; T)-{v + 2)rf*+(l 2 , l 3 , v; T) 



- {v - 2)V^ *o h, v- T) + (v - 2)y/\v-4\*t(l 2 , h, v- T)\ 

2mVE - 3 (y 3 y (v^n - (V - $+) (i 2 , i 3 , v, t) 



+ 



+ 



(4vt7\/A) 1 
16a 2 A(l + 7 2 ) 



(47T7\/A) 2 / 3 



(3.39) 



where ^ Q 4 are defined as follows: 



*±{l 2 ,h,v;T)=* 



v±A , v±2 , 

2 , l 3 ,v±A;T) + # 

v±2 1 

v±2 



v ± 4 

h,—^-h,v±4;T)+y 



v v±2 
v±2 



v±2 
v±2 



l 2 ,l 3 ,v±4;T 
■l 2 ,l 3 ,v±4;T 



h,v±4:T 



(3.40) 



and 



v±2 v 
• hi 



v±2 



v±2 



■ h,h,v,T 



v , v±2 , 

■h, h,v;T) +^ 



v±2 



v 



V 



v±2 



■ l 2 ,h,v;T 



w±2 



(3.41) 
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while ($ — $ + ) is given by 

($- - $+) (l 2 , l 3 , v- T) = (v>-2 + sgn(t;-2)| - ^ - 2 - sgn(t; - 2)|) 

x (l + sgn(v-2))*(l 2 ,l 3 ,v-2;T) 
-2(V^T3-V^TT)^(h,h,v + 2;T). (3.42) 

(The imaginary coefficients in (I3.39j) come from the action of single sin/ijCj terms.) 

Eq. (13.391) immediately implies that, as in the Bianchi I model, the steps in v are uniform: 
the argument of the wave function only involves v — 4, v — 2, v, v + 2 and v + 4. Thus, there is 
a superselection in v. For each e G [0, 2), let us introduce a lattice C e of points v = 2n + e if e 
is or 1 or v = n + e otherwise. 7 Then the quantum evolution — and the action of the Dirac 
observables Pt and V\t commonly used in LQC — preserves the subspaces T^phy consisting 
of states with support in v on C t . The most interesting lattice is the one corresponding to 
e = since it includes the classically singular points v = 0. 

Finally, it is obvious from f!3.39j) that in the limit a — > quantum dynamics of the Bianchi 
II model reduces to that of the Bianchi I model discussed in Jf]}. In particular, it is possible 
to obtain the LQC dynamics for the k=0 FRW cosmology from this model by first taking 
a —>■ and then following the projection map defined in section IVA in |6(. 



IV. EFFECTIVE EQUATIONS 

In the isotropic models, effective equations have been introduced via two different ap- 
proaches — the embedding and the truncation methods. Both start by regarding the space 
of quantum states as an infinite dimensional symplectic manifold — the quantum phase 
space — which is also equipped with a Kahler structure that descends from the Hermitian 
inner product on the Hilbert space. In the first method, one finds a judicious embedding of 
the classical phase space into the quantum phase space which is approximately preserved by 
the quantum evolution vector field 27, 28[. By projecting this vector field into the image 



of the embedding one obtains quantum corrected effective equations. In the isotropic case 
these effective equations provide an excellent approximation to the full quantum evolution 
of states which are Gaussians at late times, even in the A ^ as well as k=±l cases where 
the models are not exactly soluble. In the second method one uses expectation values, un- 
certainties, and higher moments to define a convenient system of coordinates on the infinite 
dimensional phase space. The exact quantum evolution equations are then a set of coupled 
non- linear ordinary differential equations for these coordinates. By a judicious truncation of 
this system one obtains effective equations containing quantum corrections [3l[ . In its spirit 
the first method is analogous to the 'variational principle technique' used in perturbation 
theory, in that it requires a judicious combination of art (of selecting the embedding) and 
science. It is often simpler to use and can be surprisingly accurate. The second method is 
more systematic, similar in our analogy to the standard, order by order perturbation theory. 
It is also more general in the sense that it is applicable to a wide variety of states. In this 
section we will use the first method to gain qualitative insights into leading order quantum 
effects. 



7 The lattice for e ^ 0, 1 is twice as large as that for e = or e = 1 due to the symmetry properties of the 
wave function. 
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To obtain the effective equations, without loss of generality we can restrict our attention 
to the positive octant of the classical phase space (where e = 1). Then the quantum corrected 
Hamiltonian constraint is given by the classical analogue of (13.161) : 



^+C e 

2 grav 



-icff 



0. 



(4.1) 



where 



,eff = _ P1P2P3 

' grav ~ 87rG 7 2 A£2 1 
1 



sin HiC\ sin p 2 c 2 + sin p 2 c 2 sin p 3 c 3 + sin p 3 c 3 sin p^ci 



8nGY 



= sin u,\C\ —I+7 



(4.2) 



Using the expressions (13.91) of /2^, it is easy to verify that far away from the classical sin- 
gularity — more precisely in the regime in which the (gauge fixed) spin connection and 
the extrinsic curvature are sufficiently small so that Ckfik 1 — the effective Hamiltonian 
constraint (14. ip is well-approximated by the classical one (12.241) . 

Since sin# is bounded by 1 for all 9, these equations imply that the matter density 
Pmatt = Pt/2V 2 = Pt/2PiP2P3 satisfies 



Pmatt — 



+ 



87i-f 2 AG£ 2 Pl 8vr 7 2 G 



x (l + 7 2 )x 2 



pi 



4 



(4.3) 



where we have introduced x := a^Jp2P?,/p\- The maximum of the expression in square 
brackets is attained at x = 2/(1 + 7 2 )v / A^pi, whence 



Pmatt 



< 



3 + (1 + 7 
8tt7 2 AG£ 2 



2W 



0.54p 



pi- 



(4.4) 



pi 



Thus, on the constraint surface in the phase space defined by (14.11) . the matter energy 
density is bounded by 0.54ppi. But this bound may be far from being optimal. In all 
isotropic models, the optimal bound on matter density was found to be 0.41ppi. In the 
Bianchi I model, available simulations by Vandersloot (private communication) show that 
the 'volume bounce' occurs when matter density is lower than 0.41ppi because there is also 
energy density in gravitational waves. It would be interesting to use numerical simulations 
to find out what happens for generic solutions to the Bianchi II effective equations. 

Finally, to obtain the effective equations for each variable, one simply takes its Poisson 
bracket with the effective Hamiltonian constraint. This gives the effective equations 



Pi = 7 1 ( — (sinp 2 c 2 + sinp 3 c 3 ) + ap 2 p 3 J cospiCi, 



Pi 



P2 



P3 



Pi { 
7P2 

Pj , 

7P3 



^sin pici + sin p 3 c 3 ) cosp 2 c 2 , 



^sin pici + sin p 2 c 2 ) cos p 3 c 3 , 



(4.5) 

(4.6) 
(4.7) 



21 



Cl 



1 

1 


'P2P3 


7 


Af 2 


+ 





2 




^ 3 C 3 





2 



( sin HiC\ sin /i 2 c 2 + sin HiCi sin /z 3 c 3 + sin /x 2 c 2 sin /^ 3 c 3 
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Ci cos /iiCi — sin /iiCi, 
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(4.10) 



In the "embedding approach" these effective equations provide the leading-order quantum 
corrections to the classical equations of motion Eqs. (12.271) - (12.321) . It would be very 
interesting to numerically test if the accuracy they display in the isotropic case for states 
which are Gaussians at late times carries over to the Bianchi II case. 



V. DISCUSSION 

In this paper, we analyzed the "improved" LQC dynamics of the Bianchi II model. As 
in the isotropic and Bianchi I cases, we chose the matter source to be a massless scalar field 
since it continues to serve as a viable relational time parameter in the classical as well as 
the quantum theory. It is again rather straightforward to accommodate additional matter 
fields in this framework. 

Our broad strategy is the same as that used in the Bianchi I model [6( . However, because 
Bianchi II models have anisotropies as well as spatial curvature, holonomies around closed 
curves are no longer guaranteed to be almost periodic functions of the connection. Hence, 
one cannot use them to construct the field strength operator on the LQC Hilbert space; a new 
conceptual and technical input is necessary to define the quantum Hamiltonian constraint 



operator. We overcame this difficulty by generalizing the strategy used so far |6|, [lOj, [11 
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16l . 1 171 l23l . [2J, [25| . Specifically, we used holonomies around open segments parallel to the 
fiducial triads °e°: to define a connection operator. This strategy is also inspired by methods 
introduced by Thiemann in the full theory [26[ . However, because of gauge fixing LQC does 
not enjoy the manifest diffeomorphism invariance of full LQG. As a consequence, in LQC one 
needs a principle to fix the 'length' of the open segment along which holonomy is evaluated. 
We required that the 'length' be so chosen that the field strength operator constructed 
from the resulting connection should agree with that constructed from holonomies around 
closed loops whenever the second construction is available. This guarantees that (apart from 
'tame' factor ordering ambiguities) the new procedure reduces to the one used in the LQC 
literature before. Moreover, the strategy of defining the Hamiltonian constraint through this 
connection operator can be used also in more general contexts. In particular, it enables one 
to overcome a conceptual limitation of the otherwise complete treatment of the isotropic, 
k=— 1 model given in [23|, 24j. More importantly, it extends to more general class A Bianchi 



models. A systematic treatment of the Bianchi IX model along the lines of this paper would 
be especially interesting. 

There is a second — but primarily technical — difference from the Bianchi I case: The 
Hamiltonian operator now contains inverse powers of p\. This was handled following a 
general method introduced by Thiemann to define inverse triad operators in LQG [26]. As 
usual, there is a factor ordering ambiguity. In the main discussion we used the simplest 
operator which has the same symmetries with respect to parity as its classical counterpart. 

After addressing these two issues, we obtained a well defined quantum Hamiltonian con- 
straint and showed that the singularity in Bianchi II models is resolved in the same precise 
sense as in the FRW and Bianchi I models. The Kinematical Hilbert space W§£ v can ^ e de- 
composed as = Ti^^ © 7~Lflg V where states in the singular subspace have support only 

on configurations with zero volume, while those in the regular sub-space have no support on 
these singular configurations. The Hamiltonian constraint annihilates states in Tl^g an d 
maps 7ifeg V t° itself. We also provided an explicit form of the Hamiltonian constraint which 
should be helpful in performing numerical simulations. 

Finally, we obtained effective equations using the "embedding method" introduced by 



Willis |27j and further developed by Taveras [28] in the isotropic case. There, although 



the assumptions made in the derivation fail in the deep Planck regime, the final equations 
provide a surprisingly accurate approximation to the full quantum evolution of states which 
are Gaussians at late times. This holds not only for the exactly soluble k=0, A = model 
but also for the much more complicated A ^ and k=±l models. It would be interesting 
to see if this phenomenon extends also to Bianchi II models. Furthermore, numerical so- 
lutions of these effective equations themselves may be of considerable interest because the 
simplest upper bound on matter density they lead to is higher than that in all other models 
studied so far, including Bianchi I. Numerical simulations of effective equations will answer 
several questions within this approximation. Is the upper bound optimal, i.e., do generic 
solutions to effective equations come close to saturating it? In the Bianchi I case, numerical 
simulations by Vandersloot (private communication) revealed that, unlike in the isotropic 
model, there are several distinct kinds of 'bounces.' Roughly, anytime a shear — or a Weyl 
curvature — scalar enters the Planck regime, quantum geometry repulsion comes into pay in 
a dominant manner and 'dilutes' that scalar, preventing a blow up. How do additional terms 
in the Bianchi II effective equations affect this scenario? Qualitative lessons from numerical 
simulations would be valuable in developing further intuition for various quantum geometry 
effects. 



23 



Acknowledgements: 

We would like to thank Gianluca Calcagni, Alex Corichi, Jerzy Lewandowski, Simone 
Mercuri, Tomasz Pawlowski and Hanno Sahlmann for helpful discussions. This research was 
supported in part by NSF grant PHY0854743, the George A. and Margaret M. Downsbrough 
Endowment, the Eberly research funds of Penn State, Le Fonds quebecois de la recherche 
sur la nature et les technologies and the Edward A. and Rosemary A. Mebus Fellowship. 



APPENDIX A: PARITY SYMMETRIES 

In non-gravitational physics, parity transformations are normally taken to be discrete 
diffeomorphisms Xi — > —Xi in the physical space which are isometries of the flat 3-metric 
thereon. In the phase space formulation of general relativity, we do not have a flat metric 
— or indeed, any fixed metric. Therefore these discrete symmetries are no longer meaningful 
(except in the weak field limit). However, if the dynamical variables have internal indices 
— such as the triads and connections used in LQG — we can use the fact that the internal 
space J is a vector space equipped with a flat metric to define parity operations on the 
internal indices. Associated with any unit internal vector v\ there is a parity operator 
U v which reflects the internal vectors across the 2-plane orthogonal to v l . This operation 
induces a natural action on triads e", the connections A\ and the conjugate momenta P" =: 
(1/87^7) Ef (since they are internal vectors or co- vectors). 



The triads e" are proper internal co- vectors. In previous references [6|, [32J , conventions 
were such that the spin connection rf turned out to be an internal pseudo vector. It was 
then natural to regard the Barbero-Immirzi parameter 7 to be a pseudo quantity so that 
the connection A l a has definite parity namely, it transforms as an internal pseudo-vector. 
This in turn led to the conclusion that P° is also an internal pseudo-vector (as one would 
expect because it is canonically conjugate to A l a ) [6j. While this is all self-consistent, these 
conventions lead to two undesirable consequences. First, in the classical theory, it is not pos- 
sible to reconstruct the triads e" unambiguously starting from the momenta P". Therefore, 
one cannot recover the space-time geometry unambiguously starting from the Hamiltonian 
theory. Second, the momenta P" are subject to a non-holonomic constraint which obstructs 
the passage to quantum theory a la LQG. However, if one sets conventions as in section 
UTAl then T i a) ^ ) A i a and P\ are all proper quantities and the two difficulties disappear (l9| . 
In the main text we have used this strategy. We now summarize the differences from the 
Appendix of 0] that it leads to. 

In diagonal Bianchi models, we can restrict ourselves just to three parity operations n,. 
Under their action, the canonical variables Cj,pj transform as follows: 

n i (ci,C 2 ,C 3 ) = (-Ci,C 2 ,C 3 ), ni(pi,P2,P3) = (-Pl,P2,Ps) , ( Al ) 

and the action of n 2 , n 3 is given by cyclic permutations. Thus, d and Pi are proper internal 
vectors and co- vectors. Under any of these maps n,, the symplectic structure (12.181) . the 
Hamiltonian (12.241) . and hence also the Hamiltonian vector field, are left invariant. This is 
just as one would expect because Hj are simply large gauge transformations of the theory 
under which the physical metric q a b and the extrinsic curvature K a b do not change. Also, it 
is clear from the action of (IA1I) that if one knows the dynamical trajectories on the octant 
Pi > of the phase space, then dynamical trajectories on any other octant can be obtained 
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just by applying a suitable (combination of) Ilj. Therefore, in the classical theory one can 
restrict one's attention just to the positive octant. 

Let us now turn to the quantum theory. We now have three operators Ilj. Their action 
on states is given by 

n 1 *(Zi,z 2 ,z 3 ) = *Hi,W3), (A2) 

etc. What is the induced action on operators? Since 



we have 



n^n^, z 2 , h) = it [it *(-/!, z 2 , i 3 ] 
ftififti = -k. 



(A3) 



(A4) 



The Hamiltonian constraint operator, modulo factor ordering which is not important here, 
is given by Eq. (13.161) . To calculate its transformation property under parity maps, in 
addition to (IA4I) . we also need the transformation property of the operators sin /ijCj and 
e and operators corresponding to inverse powers of p±. Due to the symmetries of type A 
Bianchi models, to know the properties of sin under parity transformations, it is sufficient 
to calculate Ilj sin/ziCillj. We have: 



nisin/!iCini^(Zi,/2,Z 3 ) = III 



2^3 | 



l 

2? 



-sinjLtiCi*(Zi,Z 2 ,Z3 



Z 2> Z 3 )-*(ii- 



|z 2 z 3 



IZ2Z3 
z 2 ,z 3 )' 



whence 



III si n Z^iCiIIi 
An identical calculation shows that 

1 



n 2 sin yUiCiITj Z 2 , Z 3 ) = -r Il 2 

2i 



sin /X1C1. 



' -Z 2 ,Z 3 )-*(Zi+ 1 



(A5) 
(A6) 



|Z 2 Z 3 | 
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|Z 2 Z 3 | 

sin/2iCi^(Zi,Z 2 ,Z 3 



j z 2 , z 3 



*(Zi 



|z 2 z 2 



|Z 2 Z 3 

") ^2) Z 3 ) 



-k,h) 



and similarly for f[ 3 . Therefore, we have: 

fl 2 sin ftiCiU.2 = sin ftici, and fl 3 sin fiiCtU^ = sin Ji\C\. 



(A7) 



(A8) 



As expected, these transformation properties of sin \X\C\ under Ilj mirror those of c\ under 
the three parity operations Ilj in the classical theory. (Note that, because of the absolute 
value signs in the expressions (13. 9p . /2j do not change under any of the parity maps.) Finally, 
it is clear from Eq. (13.81) that 



IT-Sll, 



e if v = 0, 
— e otherwise, 



(A9) 
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and from Eq. (13.251) that 

n.ipTR^HPi'R 4 . (Aio) 

(Note incidentally that this need not be the case for different factor ordering choices in Eq. 

(TOO 

We can now collect these results to study the transformation property of the Hamiltonian 
constraint. Consider first the regular subspace TCf™ v of spanned by states which have 
no support on points with v — 0. From Eq. (I3.16P it follows that the restriction to TC^ V of 

the gravitational part of the Hamiltonian constraint is left invariant under Ilj. Since f>\ is 
manifestly invariant, on the regular sub-space we have 

% C H % = Ch (All) 

Next, since the gravitational part of the Hamiltonian constraint annihilates the states in the 
singular sub-space (i.e. those with support only on those points at which v — 0), we have 

C H ^ = -h 2 dl^ = fli C H IW (A12) 

Thus, the Hamiltonian constraint operator is left invariant by all the parity operators, mir- 
roring the behavior of its classical counterpart. 

This invariance implies that, given any state ^ G , the restriction to the positive 
octant of its image under C grav determines its image everywhere on H|™ v . This property 
simplifies calculations and was used to arrive at the form of the Hamiltonian constraint given 
in fl3T39|) . 
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